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ABSTRACT 

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the 
distance between the peaks a measure of the rotation velocity. However, double-peaks can 
arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact 
solutions of the line transfer problem we present a detailed study of line emission from geo- 
metrically thin Keplerian disks. We derive the conditions for emergence of kinematic double 
peaks in optically thin and thick disks, and find that it is generally impossible to disentangle 
the effects of kinematics and Une opacity in observed double-peaked profiles. Unless supple- 
mented by additional information, a double-peaked profile alone is not a reliable indicator of a 
rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators 
of rotation in optically thick disks. 

Key words: radiative transfer — line: formation — line: profiles — methods: numerical — 
ISM: lines and bands — galaxies: active — quasars: emission lines 



1 INTRODUCTION 

Rotating disks arc key elements in numerous astrophysical sys- 
tems. Protoplanetary disks around young stellar objects (YSOs) 
are where planets eventually form, and thus are of paramount im- 
portance in all theories of planet formation. Mass transfer in cata- 
clysmic variables (CVs) proceeds through an accretion disk around 
the white dwarf. In active galactic nuclei (AGNs), the accretion 
disks are the conduit of material to the central black hole and the 
source of continuum emission at all wavelengths shorter than IR. 
Unfortunately, small angular sizes preclude direct imaging of the 
disks in most CVs and YSOs, and virtually all AGNs; only a hand- 
ful of AGNs provide direct evidence for Keplerian disks through 
interferometric observations of water masers (see Greenhill 2007). 
In the absence of imaging, a double-peaked line profile is often 
taken as evidence for an inclined rotating disk, with the distance 
between the peaks a measure of the rotation velocity. For example, 
the detection of double-peaked profiles has been used to establish 
the presence of disks in AGNs (see Eracleous et al. 2009, and refer- 
ences therein), and the analysis of double-peaked profiles has been 
used to deduce the mass of the central object from the derived Ke- 
plerian rotation (Zhang et al. (2008)). 

In the widespread reliance on the double-peaked profile as a 
marker of rotating disks it has been overlooked that such a profile 
may actually result not from disk rotation but, instead, from purely 
radiative transfer effects in optically thick, non-rotating systems. As 
has long been known, the simplest line transfer problem involving 
a uniform, quiescent slab without any large-scale motions leads to 
double-peaked profiles at large optical depths (sec Avrctt & Hum- 
mer 1965, and references therein). Unfortunately, the existence of 



these non-kinematic double peaks has been largely overlooked. In- 
stead, the Home & Marsh (1986) seminal work introduced flat- 
top line profiles for the optically thick segments of rotating disks, 
an approximation that ignores the spatial variation of the source 
function. A number of recent detailed calculations of line emission 
from protoplanetary systems did produce double-peaked profiles in 
some optically thick face-on disks (Pavlyuchenkov et al. 2007; Cer- 
nicharo et al. 2009; Ceccarelli et al. 2010), but the implications of 
these results were not discussed. The significance of non-kinematic 
double peaks has been noted only relatively recently by Hummel 
(2000) in the context of B[e] stars. 

Our goal here is to present a general study of line emission 
from rotating disks to determine whether and when the effects 
of kinematics and line opacity can be disentangled in observed 
double-peaked profiles. Following Home & Marsh we consider 
a geometrically-thin Keplerian disk and compute its line emis- 
sion profile for a large range of optical depths with detailed ra- 
diative transfer calculations. The radiative transfer problem for a 
non-rotating geometrically thin uniform disk is identical to that for 
an infinite uniform slab when neglecting radiation propagation in 
the plane. In §2 we demonstrate the emergence of non-kinematic 
double peak profiles, without any large scale motions, in the line 
emission from optically thick uniform slabs; in Appendix A wc de- 
rive some simple analytic expressions that approximate adequately 
the exact calculations. In §3 we present the solution for Keplerian 
disks, which shows double peak line emission from face-on view- 
ing. We study the properties of kinematic double peaks, which are 
introduced by the disk rotation, and compare them with the non- 
kinematic ones. Section §4 contains a discussion. 
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Figure 1. Frequency profiles 'i'(x) of fine radiation from a uniform slab, 
normalized to a unity peak flux, from exact numerical calculations with 
the CEP method. Each panel corresponds to a different value of the photon 
destruction probability e (see eq. 1). Dift'erent lines correspond to different 
values of tj, the line-center slab optical depth, as marked. 

2 NON-KINEMATIC DOUBLE PEAKS 

To demonstrate the formation of double-peaked profiles in the ab- 
sence of any large scale motions we start with the simplest prob- 
lem involving the smallest possible number of free parameters — a 
uniform slab with constant temperature and density and no radi- 
ation other than line radiation generated internally by collisions. 
We solve the problem with the Coupled Escape Probability (CEP); 
this is an exact method we developed recently for the line transfer 
problem that offers great speed advantages over traditional compu- 
tational approaches (Elitzur & Asensio Ramos 2006). We consider 
line emission from a transition between upper level u and lower 
level / separated by energy £„/ = hvg. The transition properties are 
fully specified for line transfer calculations by the parameter e, the 
photon destruction probability defined through 

1 - 6 A„, ^ 

Here C,,/ and A,,/ are, respectively, the collisional de-excitation rate 
and the spontaneous emission coefficient, T is the gas temperature, 
n the density and n^rit the transition critical density corrected with 
the Boltzmann factor (see, e.g., eq. 24 in Elitzur & Asensio Ramos 
2006). Vertical position in the slab can be specified by the coordi- 
nate T, the optical depth at line center which varies from to Tj 
between the two faces. Denote by Avu the frequency width of the 
thermal motions then the optical depth at every frequency is 

V — Vq 

t(x) = tO(x), where x = (2) 

Avd 

and where <I>(x) is a profile normalized to unity at line center. The 
problem is fully specified by just two free parameters, e and Tj, 
and the line absorption profile. In the numerical calculations here 
we take O as Gaussian, exp(-x-), but the results are applicable to 
other profiles, for example Lorentzian shape. 

Figure 1 shows the frequency profiles l'(x) of the line radia- 
tion leaving the slab vertically for different sets of the free parame- 
ters e and Tt- In optically thin slabs, 'i'(x) = <I)(x). Double-peaked 
profiles emerge when tt > 10 and become more pronounced as 
e gets smaller, i.e., as the density decreases away from the criti- 
cal density. As noted in the Introduction, double-peaked profiles in 
optically thick sources are a well known radiative transfer effect 
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Figure 2. Variation of the source function, normalized by the Planck func- 
tion of the slab temperature, with optical depth into the slab for the mod- 
els shown in figure 1 . The profiles are symmetric about the slab midplane, 
which is the endpoint of the horizontal axis in each panel. Because the axis 
is scaled with the overall optical depth, the same point in each panel indi- 
cates a different distance from the surface. Note also that S/B = I implies 
a thermalized transition (T^ = T). 

(e.g., Avrett & Hummer 1965). The explanation is rather simple. 
Consider a slab whose line center optical depth to the midplane ex- 
ceeds unity. Then the frequency x^. where ^Tj<b(x^.) = 1 defines the 
line core, the spectral region \x\ < x^ in which each half of the slab 
is optically thick. Similarly, the frequency x„ where tt<I>(x„) = 1 
defines the line wings — the entire slab is optically thin at \x\ > x„. 
For the Doppler profile, the explicit expressions for these two fre- 
quencies are 

x„ = ^JlnT^ , xj = - In 2 (3) 

At the line core, |x| < x^, each half of the slab is optically thick and 
the radiation emerges from an optical distance of only ~ 1 from 
each surface. Thus the line core intensity is roughly the value of 
the source function 5 at position t = e' (see eq. A3); that is, the 
core spectral shape reflects the spatial variation of 5, or, equiva- 
lently, the line excitation temperature since 5 = BiJ^) where B 
is the Planck function (e.g., Elitzur 1992). Because of photon trap- 
ping, Tx is maximal on the slab midplane, declining toward each 
surface. So the core intensity decreases from a peak at |x| = x^, cor- 
responding to the magnitude of on the midplane, to a minimum 
at line center, where it reflects the excitation temperature at distance 
r = 1 from the surface. At the line wings, |x| > x„, the slab is op- 
tically thin everywhere and its emission is simply l(x) = 5rTO(x), 
where 5 is the value of the source function at some intermediate 
point inside the slab (see appendix A); that is, the spectral shape of 
emission at the wings follows the line profile, ^(x) = (t>(x). The in- 
tensity now decreases as x increases because the entire slab is fully 
visible and the optical depth decreases away from line center. 

Figure 2 shows the spatial profile of S/B, where B{T) is the 
Planck function of the (constant) slab temperature, vs optical depth 
into the slab for a large range of e and Tj. This spatial variation is 
directly reflected in the properties of the emerging radiation pro- 
file (fig. 1). Once the level populations reach thermal equilibrium, 
T^ = T and S /B = 1. When Tj » 1/e, the transition is thermalized 
at slab center as is evident from the figure (see also eq. A6). Ap- 
proaching the slab surface, photon escape reduces and the source 
function is decreasing, producing the intensity dip at line center. 
When Tt increases, thermalization spreads from the slab midplane 



© 2012 RAS, MNRAS 000, 1-9 



X 

II 10" 





1 ' 1 


'1 


- 


5 X x-- 


— X ^ ^ ^ - 


: -\ 


he = 0.5 


; 

- 

- 


— 


Xe=10-' 




o — 


Oe=10"^ 




- A 


AE-IQ-'^ 

1,1 





10"^ 10"' 10° 10^ 10^ 

Figure 3. Variation of the depth of the line profile central dip with slab 
optical thickness for various values of e; the horizontal axis reflects the 
problem's scaling properties. When rx a> 1/e, the depth of the dip saturates 
at a value that varies roughly as Ve (see eq. 4). 
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Figure 4. Contours of constant line-of-sight (los) velocity V for a Keplerian 
disk between radii Rj,, and Rout ■ The observer is along the V = contour. 
In the 1st quadrant, the los velocity increases along successive contours 
clockwise from the top. 



toward its surfaces, larger portions of the source iiave - T and 
the peaks flatten out. This flattening becomes more pronounced as 
6 increases, i.e., n approaches n„it. When tt < the level pop- 
ulations do not reach thermal equilibrium even at slab center, but 
photon trapping ensures that T,^, and the source function, is always 
higher there than at the surface. 

The behavior of the central dip seen in fig. 1 is readily ex- 
plained by the spatial variation of the source function. The depth 
of the dip can be characterized by the ratio of intensities at line 
center and at the peaks. Figure 3 shows the variation of this ratio 
with optical thickness for various values of e; because of the scal- 
ing properties of the problem (see eq. A6), the horizontal axis is 
labeled with eTj. For every e, the dip reaches maximal depth when 
Tt » The dip essentially disappears at e = 0.5, correspond- 
ing to n = Rent, because the source function thermalizes throughout 
most of the slab; in this regime, 5(t) is nearly flat across the line 
core (see fig. 2). In Appendix A we discuss these issues further 
and derive semi-analytic approximations that describe adequately 
the exact numerical results; in particular, fig. Al reproduces the 
line profiles quite accurately with simple analytic expressions. The 
depth of the central dip is accurately described by eq. A7; to a good 
degree of approximation, its large-Tx limit is 



I(x = 0) 

/(x = Xc) erx » 1 



(4) 



(see eq. A8). This approximation describes reasonably well the re- 
sults shown in figure 3. 

This standard problem demonstrates that double-peaked pro- 
files need not require any large scale motion because they are gen- 
erated naturally by radiative transfer effects. We will refer to these 
as non-kinematic double peaks. In very optically thick sources, the 
emergent radiation peaks at the frequency shifts x ^ ±Xc , producing 
a peak separation of 



Ay, 



peak 



(5) 



Outside the peaks, the intensity falloff traces the shape of the 
Doppler wings. Inside the peaks, the intensity falloff toward line 
center reflects the spatial decline of the source function, i.e., the 
line excitation temperature, when the surface is approached. 



3 ROTATING DISK LINE PROFILES 
3.1 Basic Theory 

Consider a geometrically-thin disk, extending from an inner radius 
to Rm = YRi„. The intensity toward angle i from the disk nor- 
mal is Iy(R, i), where R is axial radius. Then the flux observed at 
distance D and viewing angle ; (< 90°) is L,,(;)x(cos i/D^), where 



Ui) 



,(R,i)dA 



(6) 



is the disk monochromatic luminosity in direction The observed 
flux is obtained from an integration over the circular surface area 
even though the disk image is an ellipse; the reason is that projec- 
tion on the plane of the sky distorts shapes but preserves relative 
sizes of area elements because all are scaled by the common factor 
cos i. 

The line intensity of a non-rotating disk is Iy(R, i) = 
Io(R, iyViv -vo,R, i), where lo is the brightness at line center and T 
is the spectral shape of the emerging line emission. When the line 
is optically thin, *P(v) = <I>(y). Determining *P for optically thick 
lines requires a full radiative transfer calculation; a sample of the 
results of such calculations is shown in fig. 1. In general, T varies 
with i because of the variation of optical depth with inclination; it 
can also vary with R if the vertical optical depth varies with radius. 
Switching from frequency to equivalent velocity v = c(vo - v)/yo, 
the profile becomes ^(v, i), centered on v = 0. Now set the disk 
in Keplerian rotation Vmt(R)- For viewing angle /, the line-of-sight 
(los) component of the local rotation velocity of a point at radius R 
and angle 6 from the midline is 

v-XR,6,i) = VauA-^j cos 6), where "Wout = Wiot(^out) sin ; (7) 

The line emission profile is now centered on V-(R,6,i) and the 
monochromatic luminosity of a Keplerian disk is 



L{v, i) = J Iq(R, i) T(u - V:IR, 9, i]) R 



dRdO 



(8) 



This is the basic expression for the spectral shape of disk line emis- 
sion. In the case of a uniform disk the emerging profile involves 
simply the areal integration of the intrinsic emission profile T with 
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Figure 5. Lengths of iso-velocity contours for a Keplerian disk as a function 
of line-of-siglit velocity for various values of the disk width Y = RontlRm- 



the position variation introduced by tiie Keplerian rotation. We now 
present calculations for uniform disks (^^-independent Iq and 
and various forms of Although one can expect flaring and non- 
uniformity in physical disks, the uniform thin disk can be expected 
to faithfully reproduce the essence of the radiative transfer effects 
common to all disks. 



3.2 Kinematic Double Peaks — Optically Thin Disks 

The locus of all points whose emission is centered on the same 
frequency v is the constant los-velocity curve Vz = v on the surface 
of the disk. From eq. 7, for v i^Q this curve obeys 



R 



\ V I 



(9) 



Figure 4 shows the contour plots of these curves for some repre- 
sentative velocities. In the case of a rectangular emission profile 
the emission at v is proportional to the area of a strip centered 
on the corresponding contour. When the rectangular is so nar- 
row that it can be approximated with a (5-function, the emission 
becomes roughly proportional to the contour length. The shortest 
contour is the straight line at w = 0. With v increasing, as long as 

V < Woiii the contours remain open, their curvature increases and 
with it their length. The first fully closed contour is the one for Voat, 
corresponding to the smallest los velocity on the disk midline. In- 
creasing V further corresponds to smaller radial distances on the 
midline, and the contours are now shrinking, finally reaching the 
single point at the midline intersection with the inner radius where 

V - i^in (= Vant V?). The contour lengths are calculated in closed 
form in Appendix B and plotted in figure 5 for various values of Y 
using the analytic result in eq. B2. 

The variation of contour length with los velocity captures the 
essence of kinematic double-peaked profiles in Keplerian disks: 
Voui has the longest contour and thus produces peak emission. Fig- 
ure 6 shows profiles calculated from eq. 8 with a Gaussian shape 
for the intrinsic line profile T and increasing rotation velocity. The 
sequence is best understood by considering separately the emis- 
sion from the approaching (left; fig. 4) and receding (right) halves 
of the rotating disk. With increasing rotation speed, the emission 
from each half is sliding outward in v, away from t; = 0. As long 
as VaM < Awd, the emissions from both halves largely overlap and 
there are no double peaks. Increasing further the rotation velocity. 



Figure 6. Spectral shape of the emission from Keplerian disks with Gaus- 
sian line profile 'Fq = e"'" at increasing rotation velocity. The different 
curves in each panel cover a range of disk widths, as marked. These results 
should be applicable to optically thin disks (see eq. 11). 



kinematic double peaks emerge, centered on u = ±Woui, and the 
plots begin to resemble the behavior of the length contours (fig. 5). 
In particular, the profiles for Woui = 20Awd are quite similar to the 
length curves, and the depth ratio /min//max for the central dip of 
a narrow disk (ring) with Y = I.I is nearly the same as the corre- 
sponding length ratio. As the disk width increases, both the contour 
lengths and the emission profiles show quick transitions to limit ra- 
tios, albeit at different values. From analytic expressions developed 
in Appendix B we find that 



Vi'- I Y <2 



(10) 



.4/3;r 



Y>2 



Only narrow disks (i?out < ^Rm) can produce kinematic double 
peaks with a central dip deeper than 0.4. 

The Gaussian-profile results are applicable to optically thin 
lines and can be used to assess the disk parameters that will produce 
unequivocal kinematic double peaks. Since Vant = v^aARout) sin / 
(eq. 7), figure 6 shows that double peaks require not only 
^'rot(^out) > but also sin / > AwD/Wfot(^out)- And since the los 
optical depth is tt(0 = tt/cos/, the requirement of optically thin 
emission implies cos i > Tj. Taken together, the parameter region 



Tt < 1, : > 1 



< ; < arccos Tx (II) 



A-Wd Wrot(^oul) 

will always produce kinematic double peak profiles. 
3.3 Optically Thick Disks 

Optical depth effects fall into two regimes according to the rela- 
tion between the density and the transition critical density. When 
n > /iciit (f 5^ 0.5; eq. 1) the spatial variation of the source func- 
tion is sufficiently small that the dip-to-peak contrast of the non- 
kinematic double peaks is diminished. Then the effect of increasing 
optical depth is mostly to flatten the top of the intrinsic emission 
profile and increase its width by a factor of ~ Vhrrr (see §2). In 
this regime the profile Thm = 1 - e"^' '' used in the Home & Marsh 
(1986) calculations is a reasonable approximation for the exact so- 
lution of the line transfer problem. Figure 7 presents calculations 
with this profile, demonstrating the impact of line broadening with 
increasing optical thickness of the Keplerian disk. Because of the 
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Figure 7. Spectral shape of the emission from Keplerian disks with the 
Home & Marsh profile Thm = 1 -e"'^'^'- These results should be applicable 
to sources with n ^ «crit (6 > 0.5; see §3.3). The disk radial thickness 
is Sout/^in = 10: the rotation velocity ranges from Urot(Sout) on the outer 
edge to frotC^in) = VTOfroty^out) On the inner one. Rows correspond to 
various inclinations, as marked, columns show a progression of rotation 
velocities; the apparent profile narrowing between columns arises from the 
scale change of the velocity axis. In each panel, line colors correspond to 
different values of tj, the vertical optical depth, as marked, and the emission 
is normalized to a maximum of unity. 



dependence of los optical depth on i, the viewing angle now af- 
fects not only the los velocity but also the intrinsic profile shape. 
Whereas the viewing angle effect was fully absorbed into Vout in the 
two previous figures, now it must be shown explicitly. The columns 
in figure 7 show the effect of increasing rotation velocity while 
the rows show the effect of increasing viewing angle. The top row 
shows disks viewed face-on, demonstrating the profile broadening 
with vertical optical depth' . In addition to its impact on the profile 
shape, the optical depth affects also the intensity scale — emission 
from optically thin disks is reduced by a factor of Tx/cos i from that 
of optically thick disks with the same source function. 

When n < iian (e < 0.5), optical depth effects introduce 
non-kinematic double peaks in addition to broadening the emis- 
sion spectral range (§2). This region requires full solution of the ra- 
diative transfer problem, and we utilize the emission profiles 'P(x) 
from the slab solution (see fig. 1) in the fundamental relation for 
emission from rotating disks (eq. 8). Figure 8 shows the results for 
uniform disks with 6 = 10"^ and various optical depths, rotating 
at increasing velocities and viewed from different directions; these 
results are representative of the e < 0.5 domain. As in figure 7, the 
intensity is lower in the optically thin region because of the smaller 
emitting column. Additionally, the intensity at low optical depths is 
now further reduced because of the decrease in line excitation tem- 
perature, reflecting the sub-thermal level populations. The face-on 
disks shown in the top row replicate the behavior seen in fig. 1: 
Gaussian profiles when tt < 1 , non-kinematic double peaks for the 
optically thick disks. 



' The row's three panels are identical; the apparent narrowing with rotation 
velocity simply reflects the scaling of the jc-axis with ^^((/^out). 



Figure 8. Frequency profiles, normalized to unity peak in each panel, of line 
emission from uniform, flat Keplerian disks. The intrinsic emission profile 
'F(.v) is obtained from the radiative transfer solution for a slab with photon 
destruction probability e = 10"^ and optical depth as marked. All other disk 
properties are the same as in figure 7. For clarity, the results for tt = 0. 1, 1 
and 10. marked by t, were multiplied by 15, 5 and 2, respectively. 




Figure 9. Disk rotation and optically thick line einission. Shown are 
schematic illustrations of the stationary-disk profiles when (a) n > Merit 
and {b) n < ficrit. In both cases, rotation shifts to the right (left) the emission 
from the receding (approaching) half of the disk by the amount Vqm ■ In (h), 
the leading non-kinematic peaks are labeled "outer", the trailing ones "in- 
ner". The outer peaks are always moving away from each other, the inner 
ones are first approaching and then receding as Vaut increases. 



The effects of rotation on observed profiles at slanted view- 
ing are broadly similar in both figures 7 and 8 and can be read- 
ily understood with the aid of the illustrations in figure 9 that 
sketch schematically the disk intrinsic profile in each case. Intro- 
duce Wp = Awd Vln(TT/ cos 0. Then the n > «(.rit profile can be 
described with a flat top for \v\ < (figure 9a) while the intrin- 
sic n < fZcrii profile can be taken as just two peaks with centers at 
V = ±Vp and width ~ Avo (fig. 9b). The top and bottom of each 
panel illustrate the emission from each half of a stationary disk. 
When the disk is set in rotation, the profiles of the two halves slide 
against each other in opposite directions away from the center v = 
0; the emission from the approaching half (blue profile) is shifted 
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by I'out (eq- 7) to the left, from the receding half (red) to the right by 
the same amount. Full kinematic separation of the emission from 
the two halves of the disk occurs only when > + Avd and 
the trailing edge of each profile crosses the center. Thus the crite- 
rion for emergence of kinematic double structures is now 

"""'^^""'^ sin! > Vln(TT/cosO + 1 (12) 

This replaces eq. 1 1 as the condition for kinematic double peaks 
when the disk is optically thick rather than thin. At any given tt 
(> 1), this condition can be met only for sufficiently large viewing 
angles and rotation speeds. For Tt > 10 this criterion is met in the 
right column panels of figures 7 and 8, where w,ot(^out) = 9AW|3. At 
lower rotation velocities the blue and red profiles partially overlap, 
resulting in the central peak evident in the mid-colurrms of figures 
7 and 8 in the i > 40° panels. The central peak is produced for 
rotation velocities in the range - Mj-q < f out <Vf + M}-d, namely 

■\/ln{Tj / cos i) - 1 < sini < •y/hi(TT/ cos i) + 1 (13) 

and it covers the spectral overlap region \v\ < Vom -Vj, + Avo. 

Because of its greater intrinsic complexity, the non-kinematic 
double-peaked profile of e < 0.5 sources produces some additional 
structure in rotating disks. In the schematic sketch of figure 9, the 
leading peak of the emission from each half of the disk is labeled 
"outer", the traiUng one "inner". Rotation shifts the outer peaks 
away from the center and their separation increases monotonically 
with the rotation speed. These peaks are discernible in the / > 20° 
panels of figure 8 as the outermost bumps at w = +(Wout + *^p) in 
the profiles for all Tt > 10; they are always the lowest structures, 
reflecting emission from only half of the disk. The inner peaks first 
approach each other, then overlap to produce the central peak (eq. 
13) and finally move away from each other once their trailing edges 
cross the center (eq. 12). This final regime produces the quadruple- 
peak profiles evident in the right-column panels for tt > 10 and 
i > 20°. Here the intrinsic non-kinematic double-peaked profiles 
are superposed on the underlying kinematic double-peaks with 
their centers shifted to ±t'out- The kinematic separation produced by 
rotation is centered on the dip between each pair of non-kinematic 
peaks, and not on the prominent inner peaks; the markers of the 
rotation velocity are these dips. Quad-peak profiles are seen also in 
the left column of figure 8 in the / ^ 60° panels for Tt > 100. These 
quad-peaks occur at low velocities, t^out <v^- Av-q, when the lead- 
ing edge of each inner peak has not yet reached the center. The two 
doublets making up the quad are centered on ±?;p, comprised each 
of the outer peak of emission from one half of the disk and the inner 
peak from the other half. The separation between the peaks of the 
two members of each doublet is 2tiout and the separation between 
their prominent, inner peaks is liv^ - Wom)- The inner pair of peaks 
is higher than the outer one because they get contributions from the 
two underlying non-kinematic peaks. 



4 DISCUSSION 

Double-peaked line profiles can arise from kinematic motions such 
as disk rotation, collapse or bipolar outflow. In particular, rotating 
disks produce such profile shapes for rotation velocities that obey 
the relations in eq. 11 when tt < 1 or eq. 12 when tt > 1. How- 
ever, double-peaked profiles arise also in the emission from station- 
ary optically thick sources when the density is below the critical 
density. There is no obvious way to choose between these two pos- 
sibilities from the profile shape alone since non-kinematic double 



peak profiles of face-on disks with tt > 10 (fig. 8) are quite similar 
to their kinematic counterparts in optically thin slanted disks (fig. 

6) . Because there is no clear-cut distinction between the two cases, 
a double-peaked shape must be supplemented by additional infor- 
mation before it can be taken as the signature of a rotating disk. 
Such independent input generally involves the line optical depth. 
When a transition is known to be optically thin in a given source, 
double-peaked emission can be reliably attributed to kinematic ef- 
fects and eq. 1 1 then constrains the source parameters. The depth 
of the central dip is limited in that case to Im\n/Imax ^ 0.4 in all 
but narrow disks (eq. 10). A deeper central dip implies a narrow, 
ring-like emission region with R^^t < 2/?i„. 

When the possibiUty exists that a Une is optically thick the sit- 
uation is more involved because double-peaked profiles can then 
arise from either rotation or radiative transfer or both. If it is known 
additionally that n > riait then double-peaks convincingly imply a 
kinematic origin because the intrinsic profile is nearly flat top (fig. 

7) . However, when there is reason to believe that n could be less 
than Wcrit then disk rotation can be deduced only if it can be shown 
that the peak separation exceeds Avo VlnC^x / cos i) for all reason- 
able estimates of Avj) and Tt. While this is not always possible, 
the triple and quadruple profiles seen in figure 8 can be used in- 
stead to identify optically thick rotating disks. Indeed, in certain 
circumstances these complex profile shapes could be more reliable 
indicators of rotation because they cannot be produced by radiative 
transfer alone. The emergence of these profiles is a robust result 
(see §3.3, especially fig. 9) that should not depend on the simplify- 
ing assumptions in our numerical calculations. 

To demonstrate the formation of double-peaked profiles in the 
absence of any large scale motions we presented here the simplest 
line transfer problems with the minimal number of free parame- 
ters. In particular, we considered a single two-level transition with 
line excitation purely by collisions in an environment with constant 
temperature and density. However, the results remain applicable for 
a much larger set of conditions because the central dip arises from 
the decrease of the Une excitation temperature away from the cen- 
ter toward the surface of the source. This decrease arises purely 
from radiative transfer effects in our simplified models, but it is ex- 
pected in most environments with variable physical conditions be- 
cause temperature and density generally decrease outward, ampli- 
fying the effect. A line excitation temperature that instead increases 
toward the surface is possible only under rather special circum- 
stances, such as temperature inversion for example, and our results 
would not be applicable in such cases. Similarly, our calculations 
employed only the Gaussian line profile (implicitly assuming com- 
plete redistribution) but the emergence of the central dip depends 
on the fact that the line profile decreases toward the wings, not on 
the specific shape of this decrease. The dip occurs because the emis- 
sion at the Une core follows the spectral variation of S [l/<I)(x)] (see 
§2 and eq. A3), thus the Lorentzian profile would produce the dip 
just the same, albeit with a somewhat different shape. Finally, the 
disk analysis presented here was performed for simplicity in the 
thin disk approximation. Taking account of the vertical variation of 
physical parameters requires a more elaborate calculation that adds 
more free parameters (e.g. Asensio Ramos et al. 2007) but is not 
expected to affect significantly the profile shapes found here. Our 
results should remain valid as long as the viewing angle is such that 
the disk edge contribution remains negligible and the absorption of 
the continuum does not play an important role (as, for example, 
in the case of water lines: Cemicharo et al. 2009; Ceccarelli et al. 
2010). The calculations were conducted for uniform disks, with R- 
independent line-center intensity Io{R). They are applicable to other 
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situations because a disk whose /o(i?) declines steeply with R can be 
approximated with a uniform disk whose boundaries are limited to 
the peak emission region. This explains the Home & Marsh (1986) 
modeling results for the ecUpsing dwarf nova Z Cha. The double 
peaked profile in this source has /min/^max = 0.4, roughly the deep- 
est dip possible for optically thin Keplerian disks with R„^i > 2i?i„. 
While the Horne & Marsh model disk had a much larger R„„i ~ 
10i?m, it also had Iq(R) oc 1/R^. This steep intensity decline ensures 
an effectively ring-like emission region that enabled reproduction 
of the observed contrast. 

Our results show that detection claims of disks based on 
double-peaked profiles of molecular lines must be always scruti- 
nized carefully for optical depth effects. Double-peaked profiles of 
presumably optically thin lines (e.g. from rare isotopes) can be used 
reliably for identifying candidate disks around YSO whose average 
sizes are too small to be resolved even with ALMA. The case of 
AGN, where double peaked Ha and Mg II AH9S profiles were in- 
terpreted as rotating disk emission (Eracleous et al. 2009), is more 
complicated. Because their excitations are driven by photoioniza- 
tion and not collisions, our model calculations are not directly ap- 
plicable to these lines. Still, their large optical depths, ~ lO'-lO* 
(H. Netzer, private communication), imply that the source func- 
tion is unlikely to be constant in either case, triggering the non- 
kinematic double peaked shapes seen in §2. An actual estimate of 
this effect will require detailed photoionization model calculations 
coupled with exact radiative transfer calculations for these lines in- 
stead of the common escape probability approach. 



APPENDIX A: SLAB ANALYTIC APPROXIMATIONS 

Neglecting radial coupling, the radiative transfer problem of a face- 
on flat disk is identical to that of a slab, as noted above (§3). Here 
we present some approximate analytic expressions that reproduce 
adequately the solution of the slab problem. 

The two-level problem in a flat uniform slab is fully defined 
by the photon destruction probability e and the total optical depth 
Tj across the slab. Positions in the slab are characterized by the 
vertical distance t from one face, < t < tt. The formal solution 
for the intensity emerging perpendicular to the slab face is 



I(x) = <b(x) f ^ 5(T)e-<"T-^)*W^.r 
Jo 

Whenever tt<1>(x) < 1, 



Jo 



S(T)dT = STT^{x) 



(Al) 



(A2) 



where 5 is the value of the source function at some (a-priori un- 
known) point inside the slab. We approximate the emission in the 
optically thin wings (tt$(x) < 1, \x\ > x^; see eq. 3) with this ex- 
pression. In the optically thick core (Tj<t>(x) > 2, \x\ < x^), accord- 
ing to the Eddington-Barbier relation the radiation emerges from an 
optical depth of ~ 1 from either surface, therefore the intensity is 
approximately S[t = l/(b{x)]. The frequency range < \x\ < x^, 
where 2 > ttC>(x) > 1, does not lend itself naturally to any physi- 
cal approximation since the optical path is thick to one surface and 
thin to the other (except for the midplane where both paths are op- 
tically thin). However, this is a small spectral region whose relative 
thickness (x,., - x^)lxc decreases with overall optical depth (eq. 3), 
therefore we approximate the intensity there with a constant. Join- 
ing smoothly the core and wing approximations with this constant 
yields S = S(tt/2) — the unknown S is the source function on the 



slab midplane. Thus our approximation for the emergent radiation 



I{x)^ 



5[l/0(x)] 

5(tt/2) 

5(tt/2)tt*(x) 



Xc < 



\x\ ^ Xyj 

\x\ > x„ 



(A3) 



Completion of the intensity expression requires an approximation 
for the T- variation of S inside the slab. Motivated by Eq. 2.18 of 
Avrett & Hummer (1965), we fitted the source function throughout 
the slab with the following empirical expression containing the free 
parameters a, b and a: 



S{t; e,TT) = a + bkiie^T), 



(A4) 



where 



ki(r) : 



and 



K2{T) ■■ 



1 - 2 [^2(^) + KliTT ■ 



T)] 



1 

•0r 



<b{x)E2 [Ttbix)] dx. 



with E2 the 2nd exponential integral. We find that a = 0.7 yields 
satisfactory agreement with the exact numerical solutions for all 

values of 6 and tj. The constants a and b are related to the values 
of the source function on the slab surfaces and its midplane via 



S(0)-bki(0) 
5(tt/2)-5(0) 



(A5) 



The values of S (O) and 5(77/2) can be approximated utilizing some 
general properties of the slab radiative transfer problem derived in 
Avrett & Hummer (1965). They show that the source function on 
the surface of a slab with an overall optical depth tt can be related 
to that at depth tt into a semi-infinite atmosphere. For the latter 
problem there exist some general theorems in the limits of both 
large and small e. Utilizing those, we devised the following simple 
relations: 

e -I- er-T 



(A6) 



5(0) = v^V!±A(£Ii)!llMli)!!5, 

1 +j8i(6Tt)''1 +M<^tt)''^ 

where B is the Planck function of the (constant) slab temperature. 
We find that 5(0) is reproduced satisfactorily with /3i = 1.41, 
^2 = 1-23, pi = 0.6 and p2 = 0.55. This completes our analjftic 
approximation for 5(t) in all slabs. Inserting this expression in the 
approximation for the emerging intensity (eq. A3) produces the line 
profiles shown in figure Al. These profiles are calculated for the 
same parameters used in the exact model calculations presented in 
figure 1. Comparison of the two verifies the adequacy of the ana- 
lytic approximation provided by the combination of equations A3 
through A6. 

The depth of the central dip in double-peaked profiles can be 
characterized by the ratio of intensities at Une center and at the 
peaks. From eq. A3, 

/„i„ /(x = 0) 5(T=1) 



(A7) 

/max l(X = Xc) 5 (T = Tt/2) 

This expression provides an excellent approximation for the exact 
results shown in figure 3 for the variation of the central dip with 
optical depth for various values of e, but 5(t = 1) is not easily dis- 
cernible from the somewhat involved approximation for the source 
function. Neglecting the variation of 5 close to the slab surface and 
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Figure Al. Line radiation emergent from slabs with the parameters of the 
models displayed in figure 1 . These profiles were calculated analytically by 
combining the approximate expressions in equations A3 through A6. 



replacing 5(1) with 5 (0), which is given explicitly in eq. A6, yields 
the simpler, though somewhat less accurate, result 

/Oc = 0) ^ 5(r^0) ^^^^ 
I{x = x^) S(t = Tj/2) 

Combining this relation with eq. A6 produces the approximation 
given in eq. 4. 



APPENDIX B: KEPLERIAN CONTOURS 

All points on the disk surface whose emission is centered on the 
same frequency v lie on the contour defined in eq. 9. The polar 
angle varies on the contour between the inner and outer radii in the 
range [Bi^, d^^t], where 

r V 



fir 



(Bl) 



V t;oL,t < w < I'm 

Straightforward integration produces the contour length 

m = Ro.. (^)' [mj - mud] (B2) 

where 

fid) = sineVl +3 sin 6*2 + — Inf V3sinf9+ Vl + 3 sin 02 ) 
V3 ^ ' 

This relation produces the profiles shown in figure 5. Comparison 
with fig. 6 shows good agreement with the results of full calcula- 
tions with narrow-width Gaussian shape for the intrinsic line emis- 
sion. In particular, the depth ImiJImia of the central dip is nearly the 
same for the narrow disk (ring) with Y = 1.1 as the length ratio 
^min/^max ^01 the Same Y. Therefore, we can use the contour length 
for an estimate of the central dip depth contrast in narrow disks: 

/ ■ / 

'mm '^m 



'mm •'mm 
^max ^niax 



The longest contour, v ■ 
that 



, has ( 



: and ( 



(B3) 
arccos Y '-so 



Cax = i?ou,/(Varccosy-i/2) ^-^^ 2R,u, VF^ 

while the length of the shortest contour, w = 0, is simply 



(B4) 



(B5) 

Combining the last three relations produces the narrow-disk result 
listed in eq. 10. 

Although the length calculation reproduces reasonably well 
the profile shape and F-dependence also for wide disks (K » 1; 
see figures 5 and 6), it misses on the actual value of /max/^iin> the 
dip contrast; the length calculation produces accurately the central 
contrast for narrow disks but not for the thicker ones because of a 
fundamental issue with the approach to a ^-function limit. The line 
profile can be considered narrow when Avo « v, and this condition 
is never applicable at line center, where w = 0. Since the line wings 
of the closed contours contribute more to the w > contours of 
wide disks, the length calculation misses the proper result; instead, 
a calculation of the emitting area at velocity v is necessary. The 
area under the u-contour in the 1st quadrant is 

MV) = (^)' kWn) - g(0ou,)] (B6) 

where 

3 3 1 . 

s(9) = -6+ - sin6'cos(9+ - sin (9 cos 
* 8 8 4 

At large Y we have ^ njl for the Wout-contour (eq. Bl) so that 
g(^in) =7!^ ^rid the area under contours close to the longest one is 

37r 



A(v)- 



32 



r,2 / ^OUl \ 



(B7) 



For a rectangular velocity profile with width Auq <k Wom, it is 
straightforward to show from this expression that the area of peak 
emission in the 1st quadrant is 

37r 



■^peak 



: A(i;o„, - Aud) - A{v„ud - 



Ri 



(B8) 



At line center {v = 0) the emitting area in the first quadrant is be- 
tween the vertical axis and a roughly radial ray with v = Avo- This 
ray has inclination Aud/Uou, to the y-axis so that the area of dip 
emission is 



_ 1 2 

Adip — - Rai 



AVn 



Therefore, the limiting depth is 



/m 



Apeak StT 



(B9) 



(BIO) 



This result, listed in eq. 10, reproduces fairly accurately the depth 
contrast of the narrowest Gaussian profile (Avj^/Vout = 1/20) in 
figure 6. 
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